
DESCRIPTION 

METHOD FOR ESTIMATING NOISE POWER AND NOISE POWER ESTIMATION APPARATUS 
Technical Field 

5 The present invention relates to a method for estimating noise 

power and a noise power estimation apparatus, and is preferable for 
being applied to, for example, radio receiver apparatus. 

Background Art 

10 In multi carrier-code division multiple access (MC-CDMA) , there 

is a method for spreading to a frequency domain. In this method, after 
the information data stream is spreaded with a given spreading code 
series, each chip is mapped to different subcarriers. Fig. 1 is a 
schematic diagram showing a packet configuration of the MC-CDMA 

15 conducting a frequency domain spreading . Here, a primary-pilot channel 
(P-PICH) is shown with oblique lines, a dedicated-pilot channel (D-PICH) 
with halftone dot meshing, and dedicated physical channel (DPCH) with 
white background. The P-PICH is a pilot for a cell searching and a 
synchronized tracking, and the D-PICH is a pilot employed for SIR 

20 estimation for adaptive radio parameter control, channel estimation 
of the DPCH and Doppler frequency estimation. 

In these packets shown in this figure, one frame (54 symbols) 
is composed of two slots. In one frame, P-PICH, D-PICH and DPCH are 
time division-arranged. As can be seen from the figure, four symbols 

25 are arranged in one frame for the D-PICH. 

Conventionally, in order to ideally conduct a control such as 
adaptive modulation or a resource allocation, it is necessary to 
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appropriately measure the quality of received signal, and noise power 
maybe employed as a quality of received signal . The noise power includes 
a thermal noise generated inside a receiver and a noise generated by 
an interference of the other cell through the propagation path. The 
5 D-PICH having the above-described packet configuration is employed for 
the estimation of the noise power (referred to as "N+I other") per one 
frame. Conventional estimation for the noise power ("N+I other") is 
described as follows. 

Fig . 2 is a block diagram showing a configuration of a conventional 

10 receiving system. In this figure, radio processing section 11 receives 
signal transmitted from a communication partner through an antenna, 
and predetermined radio processing such as down conversion and AD 
conversion is conducted on the received signal, and radio processed 
signal is output to GI removal section 12 . The GI removal section 12 

15 eliminates guard interval of signal that is output from the radio 
processing section 11, and the processed signal is output to a FFT section 
13. The FFT section 13 acquires signal, which is transmitted on each 
subcarrier by conducting fast Fourier transform for signal output by 
the GI removal section 12 . The acquired signals at each of subcarriers 

20 are output to a pilot extract section 14. 

Pilot extract section 14 extracts time multiplexed pilot channel 
(PICH) from a signal of each of subcarriers that is output from the 
FFT section 13, and the extracted PICH is output to a noise power 
estimation section 15. On the other hand, signal excluding PICH is 

25 output to P/S converter 16. 

The noise power estimation section 15 estimates noise power on 
the basis of PICH that is output from the pilot extract section 14. 
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Details of the noise power estimation section 15 will be described later . 

P/S converter 16 conducts P/S conversion processing for signal 
excluding PICH that is output from the pilot extract section 14 and 
output P/S converted signal toadespreading section 17 . Thedespreading 
5 section 17 conducts a despread processing by multiplying a predetermined 
spread code to the P/S converted signal and output despreaded signal 
to a demodulating section 18. The demodulating section 18 conducts 
demodulating processing for despreaded signal that is output from the 
despread section 17, and output them to a decoding section 19. The 

10 decoding section 19 conducts a decoding processing such as turbo decoding 
and the like to signal that is output from the demodulating section 
1 8 , and output them to an error detecting section 2 0 - The error detecting 
section 20 conducts an error detection of signal that is output from 
the decoding section 19, and output data if error is not detected. 

15 Fig. 3 is a block diagram showing an internal configuration of 

the noise power estimation section 15 in a conventional apparatus . Adder 
21 calculates pilot symbol correlated/in-phase added value (^i) of a 
subcarrier i in time orientation using PICH that is output from the 
pilot extract section 14, and output the calculation results of each 

20 of subcarriers to a delay device 22-1 (shown as "D" in figure) and a 
multiplier 23-1. 

The delay device 22-1 delays ^ of respective subcarriers that 
are output from the adder 21 until next ^ is input, and the obtained 
^ is output to a delay device 22-2 and amultiplier 2 3-2 . These operations 

25 are similarly applied to delay devices 22-2 to 22-4, and these delay 
devices delays to ^i+i that are output from the adder 21 (i is defined 
to be within a range of from 0 to 7 67) by one, respectively, to output 
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H the delay devices . However, there 
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orientation (5) is obtained for respective subcarriers by using ri,j. 
More specifically, the values are obtained by the following formula: 
[formula 1] 




5 Alternatively, E,± can be represented by = hi + ni, by using 

noise component ni and channel variation component hi. ^i calculated 
in the adder 21 is transmitted to delay devices 22-1 and multiplier 
23-1 consecutively. 

We will specifically find out "N+I other" by focusing on subcarrier 
10 i as follows. Here, a case of using ^ averaged for two adjacent 
subcarriers for calculating noise power per one subcarrier will be 
described: 
[formula 2] 

15 [formula 3] 

Here, Ci-o.s and Ei+o.s are the averaged values of respective ^ of 
subcarrier number i-1 and i and subcarrier nxjmber i and i+1. Assuming 
that the condition is free of noise and the correlation between adjacent 
20 subcarriers is 1, the value obtained by subtracting Ci+o.s from Ci-o.s 
becomes zero. Therefore it is considered that the values obtained as 
differences are noise component. The noise power for subcarrier i is 
found by squaring the noise component thus obtained . More specifically, 
it can be expressed by the following formula (4) : 



[formula 4] 

The following formula (5) can be derived from formula (2) and 
formula (3) : 
[formula 5] 



In this case, correlations between adjacent subcarrier waves 
are not necessarily be 1 for subcarriers i-1 to i+1, and thus there 
may be a case that a difference in channel variation h is generated. 
In other words, since the value obtained by subtracting hi+i from hi_i 
is not zero, error is generated for ai that is wanted to be found. 
Therefore, a correction is done in order to diminish error generated 
by channel variation to obtain the following formula (6) : 
[formula 6] 



Here, h'^i-i is a value obtained by averaging \ for the plurality 
of subcarriers adjacent to subcarrier i-1, and can be obtained by the 
following formula (7) : 
[formula 7] 
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(5 




(6 





= -l^i-i + 2/j,., + A, + «(-2 + 2«,_, + «, ) • • • (7 



On the other hand, h'^i+i is obtained by the following formula 

(8) : 

[formula 8] 



^i+l ~ ^ (^1+0.5 "^^1+1.5) 



2 

= + 2/1,^, + + '^z + 2w,^, + ) ... (8 
4 

In this way, h can more precisely be found by conducting the 
averaging operation using ^ across a plurality of subcarriers, thereby 
diminishing an error of ai. Formula (7) and formula (8) are substituted 
into formula (6) to obtain formula (9) shown below: 
[formula 9] 



= 



Here, assuming that differences of respective channel variation 
components are equal, i.e., hi_i - hi+i = hi-2 - hi = hi - hi+2 = A, the 
relationship of the following formula (10) is satisfied: 

[formula 10] 

2 A = 2(/i,.. - /I,,, ) = {h,_, -/!,) + {h, - A,,, ) ... (10 

Formula (11) is provided by conducting an approximating 
calculation of formula (9) using formula (10) : 
[formula 11] 

Formula (11) presents an arithmetic operation conducted in the 
square device 25, and, in practice, calculation is carried out by 
replacing n with ^ . 
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[formula 12] 

a/ =[|J|-^« +2#,., -2^,„ -.. (12 

In this formula, coefficients of respective E, are the coefficients 
that are set by multipliers 23-1 to 23-4, respectively. In other words, 
5 a part of the calculation using formula (12) is conducted in these 
multipliers , The results by the calculation within the square device 
25 is output to the cumulative adder 26. Thus, provided that ni-2 to 
ni+2 satisfy probability distribution of Gaussian distribution in formula 
(11) and average power is presented as (na)^, it can be expressed by 
10 formula (13) : 
[formula 13] 

^;=(^J-9(«j^ - (13 

Assuming that a = (1/8)^ (9), the noise power ("N+I other") can 
be presented as: 
15 [formula 14] 

In this case, cumulative-adding operation for a' that is output 
from the square device 25 is conducted in the cumulative adder 26, and 
4/ ( (Nc -4) X a) are multiplied thereto in the multiplier 27, thereby 
20 providing an estimation of "N+I other" - Here, the reason of multiplying 
4 thereto is to reflect the value into noise power per one frame ("N+I 
other") / as ^ is obtained by conducting the averaging operation using 
four pilot symbols. In addition, the reason of dividing by (Nc-4) is 
that number of cumulati ve -added a' is given by {Nc-4) , and thus conducting 
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the average operation is intended. 

Estimations of noise power has conventionally been carried out 
by the method stated above. 

However, in the above-described conventional process for 
5 estimating noise power, there is a problem, in which noise power cannot 
be precisely estimated when the channel variations in respective 
subcarriers are large by the influence of frequency selectivity fading 
or the like. 

10 Disclosure of Invention 

An object of the present invention is to provide a method for 
estimating noise power and a noise power estimation apparatus, which 
provide improvement in the estimation accuracy of the noise power, even 
if the correlation between adjacent subcarriers is low due to the 
frequency selective fading. 

The above-described object is achieved by repeatedly correcting 
differences in the channel variations in different subcarriers for 
multiple times when noise power per one frame ("N+I ©ther") is estimated. 

Brief Description of Drawing 

Fig . 1 is a schematic diagram of packet configuration of MC-CDMA; 
Fig . 2 is a block diagram showing a configuration of a conventional 
25 receiving system; 

Fig. 3 is a block diagram showing an internal configuration of 
a conventional noise power estimation section; 
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Fig, 4 is a diagram showing PICH with its relation between time 
and subcarrier; 

Fig. 5 is a block diagram showing an internal configuration of 

a noise power estimation section according to the embodiment 1 of the 
5 present invention; 

Fig. 6 is a block diagram showing a configuration of receiving 

system according to the embodiment 2 of the present inventions- 
Fig . 7 is a block diagram showing an internal configuration of 

a noise power estimation section according to the embodiment 2 of the 
10 present inventions- 
Fig. 8 is a block diagram showing an internal configuration of 

a noise power estimation section according to the embodiment 3 of the 

present invention; 

Fig. 9 is a block diagram showing an internal configuration of 
15 a noise power estimation section according to the embodiment 4 of the 

present invention; and 

Fig. 10 is a block diagram showing configuration of receiving 

system according to the embodiment 5 of the present invention. 

20 Best Mode for Carrying Out the invention 

Embodiments of the present invention will be described in 
reference with the annexed figures as follows. 

25 (First Embodiment) 

Fig. 5 is a block diagram, showing an internal configuration 
of a noise power estimation section according to the first embodiment 
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of the present invention. However, in this figure, same numeral number 
as in Fig. 3 is assigned to an element that is commonly appeared in 
Fig. 3, and the detailed explanation thereof is omitted. The point 
that Fig. 5 is different from Fig. 3, is that the number of the delay 
5 devices and themultipliers are respectively increased, and coefficients 
for multipliers 202-1 to 202-6 and a multiplier 203 are changed. 

An adder 2 1 calculates pilot symbol correlated / in-phase addition 
value (^i) in time orientation of subcarrier i using PICH that is output 
froma pilot extract section 14 , and the calculated results for respective 
10 subcarriers are output to a delay device (indicated by "D" in the figure) 

201- 1 and a multiplier 202-1. 

Delay devices 201-1 to 201-6 delays ^1-3 to ^1+2 (i = 0 to 767) 
that are output from the adder 21 by one, respectively, and the delayed 
5 is output to the multipliers that are respectively coupled to delay 
15 devices. However, there is no output from the delay device 201-3 to 
the multiplier. 

Predetermined coefficients are set in advance, and the multipliers 

202- 1 to 202-6 multiply predetermined coefficients to ^ that are 
respectively input, and the multiplied values are output to the adder 

20 24. 

The multiplier 203 multiplies 4/((Nc - 6) a) to values that are 
output from a cumulative adder 2 6 to calculate "N+I other"* Here, Nc 
is the number of the employed subcarriers, and a is a constant. 

Next, an arithmetic operation conducted in a noise power 
25 estimation section stated above will be described using formula. 
However, since the processes of foinnula (1) to (9) are same as stated 
above, description will be made for operations after formula (9) . In 
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formula (9), a correction same as formula (6) is further introduced 
in order to diminish an error caused by channel variation component 
h. This is because, when a drift between the subcarriers is occurred 
due to frequency selective fading, correlation of the adjacent 
5 subcarriers is decreased to increase values obtained by subtracting 
the channel variation components generated between subcarriers 
therefrom, and thus the purpose for the correction is to estimate noise 
power ("N+I other") with higher accuracy. The introduction of the 
correction into formula (9) provides a formula (15) shown below: 
10 [formula 15] 



-«/-2 +2«,-i -2«,^, +«,^2r 

Expand formula (15) to obtain formula (16) : 
[formula 16] 



15 C7, = 



Here, assuming that differences of respective channel variation 
components are equal, i.e., hi-i - hi+i = hi-3 - hi-i = hi-2 - hi = hi - 
hi+2 = hi+i - hi+3 = A, the relationship of the following formula (17) 
20 is satisfied: 
[formula 17] 

6A = 6(V, -/»,.,) = -(V3 - V.)+4(/,,-, -/i,)+4(/i, -V2)-(A,.. -^,.3) - (17 

Approximating formula (16) by employing formula (17) provides 
the following formula (18) : 
25 [formula 18] 
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Actually, the calculation of this formula (18) is carried out 
by replacing n with ^. 
[formula 19] 

v2 



r" |#,-3 -4^/-2 +5^,,, -5^,,, +4^,,, -^,,3p 



(19 



Coefficients set for multipliers 202-1 to 202-6 are same as 
coefficients of respective ^ presented by formula (19), and values of 
the coefficients reflect the twice-corrections - Further, the adding 
operation for ^ is conducted in the adder 24, and the added value is 
10 squared in a square device 25. Here, provided that ni_3 to ni+3 satisfy 
probability distribution of Gaussian distribution in formula (18) and 
the average power is presented as (na)^, it can be expressed by formula 
(20) : 

[formula 20] 



^l-]\fi4{nj ... (20 



U2> 

Assxaming a = (1/32)^ (84), noise power ("N+I other") is presented 

as : 

[formula 21] 



hiker =-\<^3 +^4 +^5 ••' + ^A^.-4 J/^c-"6 



(21 



20 In the cumulative adder 26, cumulative-adding operation is 

conducted for a" that is output from the square device 25, and 4/ ( (Nc 
- 6) a) are multiplied thereto in the multiplier 203. (Nc -6) is number 
of the cumulative-added a" . 
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As such, according to the present embodiment, errors generated 
by the channel variations can be diminished by introducing corrections 
twice during the process for estimating noise power, and thus the 
estimation accuracy of the noise power can be achieved, even if the 
5 correlations between adjacent subcarriers are decreased due to the 
frequency selective fading. 

(Second Embodiment) 

While the first embodiment illustrates the case of using ^ that 
10 are obtained by averaging ^ of adjacent two subcarriers when noise power 
per one subcarrier is calculated, the second embodiment will illustrate 
a case of using ^ that are obtained by averaging ^ of adjacent four 
subcarriers . 

Fig. 6 is a block diagram showing configuration of receiving 
15 system according to the second embodiment of the present invention. 
However, same numeral number as in Fig. 2 is assigned to an element 
that is commonly appeared in Fig . 2 , and the detailed explanation thereof 
is omitted. A delay spread determination section 201 generates a delay 
profile on the basis of signal (received signal) that is output from 
20 a radio processing section 11 to measure a level of a spreading of 
multipath (Delay Spread) . 

The measurement results represent conditions of frequency 
selective fading, and levels of the correlations of the adjacent 
subcarriers can be judged by using thereof. 
25 More specifically, the delay spread determination section 201 

determines levels of the correlations of the adjacent subcarriers based 
on the measurement results, and output the evaluation results to the 
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noise power estimation section 202 . 

The noise power estimation section 202 changes number of ^ employed 
for the averaging operation, based on the evaluation results by the 
delay spread determination section 201. In other words, when the 
5 correlations between the subcarriers are high, the accuracy for the 
estimation of the noise power is improved by increasing number of ^ 
employed for the averaging operation. On the other hand, when the 
correlations between subcarriers are low, an increase of the error caused 
by the differences of the channel variations can be prevented by reducing 

10 number of ^ employed for the averaging operation. In the second 
embodiment, as compared with the first embodiment, the correlations 
between subcarriers are set to be higher, and number of ^ employed for 
the averaging operation is increased to 4. 

Fig. 7 is a block diagram showing an internal configuration of 

15 the noise power estimation section 202 according to the second embodiment 
of the present invention. In this figure, same numeral number as in 
Fig. 5 is assigned to an element that is commonly appeared in Fig. 5, 
and the detailed explanation thereof is omitted. The point that Fig. 
7 is different from Fig. 5, is that the adder 21 is replaced by an adder 

20 301, number of the delay devices and the multipliers are respectively 
increased, and coefficients for multipliers 303-1 to 303-12 and a 
multiplier 304 are changed. 

Arithmetic operations provided by the noise power estimation 
section 202 shown in Fig. 7 will be described using formula. Here, 

25 specific arithmetic operations will be described focusing on a 
subcarrier i: 
[formula 22] 
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^/-o.5 =\i^i-2 +«/-2)+;J;(^-i +'»i)+;J;(^,>i 



(22 



[formula 23] 



4 4 4 4 

Similarly as in the first embodiment, a difference value obtained 
by subtracting €±+0.5 from ei-o.s is squared to obtain noise power for 
the subcarrier i. 
[formula 24] 

^1 =K-o.5 - -^1+0.5 r (24 

Formula (22) and formula (23) are substituted into foarmula (24) 

and the resultant formula is further simplified to obtain the following 

formula (25) : 

[formula 25] 

2 



7(^i-2 -'^.>2)+7(«i-2 -«i>2) 

4 4 



(25 



Here, a correction is added to formula (25) in order to diminish 
error generated by channel variation: 
[formula 2 6] 

2 



= 



^{^-2 -^,>2)+;^(^i-2 -«i>2)-;^(4-2 -^^i) 



(26 



h'^i-2 and h''i+2 are presented by the following formula (27) and 
formula (28), respectively: 
[formula 27] 

4-2 =^(^/-2.5+^i>1.5) 
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= ^('',-4 + 2/»,.3 + 2h,.2 + 2/»,_, + h, + n,_, + 2n,.3 + 2/i,_, + 2n,., +«,)••• (27 
[formula 28] 

* 1 / \ 

Formula (27) and formula (28) are substituted into formula (26) 
to obtain formula (2 9) : 
[formula 29] 

10 = j^±J|{-(/^,,^, ^/,J_2(V3 -/i,.^,)+6(V, -/^,>2)-2(/i,_, -A..3)-(/^. -/i,>4)} 

+ {- n,_, - 2n,.3 + 6a7,_2 - 2«,_, + 2/7,,, - 6/7,,2 + 2«,,3 + «, }|' • • • (29 

A correction is further introduced to formula (2 9), for the purpose 
of reducing an influence thereto by the channel variation components, 
in consideration of frequency selective fading. 
15 [formula 30] 

<y- = (^J|{-('^'-4 -/i.)-2(/i,.3 -^..,)+6(/i,., -/i,,J-2(v, -K,y{h-h,,,)\ 

- {-(V4 -4)- 2(^3 -4,J+6(4., -4.2)-2(v, -K^V% 
+ {- «,_4 - 2w,_3 + 6«,,2 - 2«,-_, + 2n,,, - 6w,„2 + 2«,,3 + n,^^ }|' • • • (30 
Expand formula (30) to obtain the following formula (31) : 
20 [formula 31] 
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-20(/,,., -/»,,3)-8(/», -/»,.J+4(/,,„ -h,J+ih,,, -h,J} 

+ {n,^ + 4«,._j - 8/i,._4 - 20«,_3 + 45/i,_2 - 24«,_, + 24n,^, - 45n,^2 

+ 20/i,,3+8«,^-4»,,,-«,^}|' - (31 

Here, assuming that differences in respective channel variation 
components are equal, the following formula (32) is obtained: 
[formula 32] 

46A = 46{h,., -h,,,)=-{h,_, - V,)-4(/7,.3 -h,.,)+S{h,^ -h,)+20{h,_, -h,,,) 
+ 20{h,_,-h,,,)+S{h, -h,J-4ih,,, -h,J-{K2 -h.J - (32 

Approximate formula (31) by employing formula (32) to obtain 
the following formula (33) : 
[formula 33] 

" f 1 VnV 

o-/ = I —J g J h/-6 + 4«,.5 - 8«/-4 - 20«,_3 + 45«,_2 - 24/1,,, + 24«,^, - 45n,^^ 

+ 20«,,3+8«,-^-4m,-,3-«,^|' ... (33 

Actually, the calculation of this formula (33) is carried out 
by replacing n with ^ . Calculating operation for obtaining square 
thereof is conducted in the square device 25: 
[formula 34] 

v2 



+ 4^,.3 -8^,^ -20^,.3 +45^,., -24^,., +24^,,, -45^,, 



+ 20^,,3+^,,,-4^,,3-^,,,|^ - (34 

Coefficients set for multipliers 303-1 to 303-12 are same as 
coefficients of respective ^ presented by formula (34), and values of 
the coefficients reflect the twice corrections. Further, adding 
operation for ^ is conducted by the adder 24, and the added values are 
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squared in the scjuare device 25. 

Here, provided that ni-6 toni+6 satisfy probability distribution 
of Gaussian distribution in formula (33) and average power is presented 
as (na)^/ it can be expressed by formula (35): 
[formula 35] 

i2 



^_L) .a,64„.' ... (36 

Assuming a = (1/256)^ (6164), noise power ("N+I other") is presented 



as : 

[formula 36] 



4 / " " " / 



(36 



In the cumulative adder 2 6, cumulative-adding operation is 
conducted for o" that is output from the square device 25, and 4/ ( (Nc 
- 8) a) are multiplied thereto in the multiplier 304. (Nc - 8) is nxamber 
of the cimiulative-added a" . 

15 As such, according to the present embodiment, errors generated 

by the channel variations can be diminished by introducing corrections 
twice during the process for estimating noise power, when ^ of adjacent 
four subcarriers are averaged in the calculation for obtaining noise 
power per one subcarrier, and thus the estimation accuracy of the noise 

20 power can be achieved, even if the correlation between adjacent 
subcarriers is decreased due to the frequency selective fading. 



(Third Embodiment) 

While the first embodiment illustrates the case of twice 
25 conducting the correcting operations during the process for estimating 
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noise power, the third embodiment will illustrate a case of conducting 
the correcting operations for three times during the process for 
estimating noise power. 

Fig. 8 is block diagram showing an internal configuration of 
5 a noise power estimation section according to the third embodiment of 
the present invention. In this figure, same numeral number as in Fig. 
5 is assigned to an element that is commonly appeared in Fig. 5, and 
the detailed explanation thereof is omitted. The point that Fig. 8 
is different from Fig. 5, is that the number of the delay devices and 

10 the multipliers are respectively increased, and coefficients for 
multipliers 402-1 to 402-8 and a multiplier 403 are changed. 

Arithmetic operations provided by the noise power estimation 
section shown in Fig . 8 will be described using formula . Here, specific 
arithmetic operations will be described focusing on a subcarrier i. 

15 However, since the operations until formula (16) is derived in this 
example is the same as in the first embodiment, description will be 
made for operations after obtaining formula (16) . 

In the third embodiment, additional correction is further added 
to formula (16), and thereafter, similarly as each of the embodiments 

20 stated above, expansion and approximation are carried out to obtain 
the following formula (37) : 
[formula 37] 

Actually, the calculation of this formula (37) is carried out 
25 by replacing n with ^. 
[formula 38] 
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<=[:^]'|-^,-4 +6^,.3 -14^« +14^,-, -14^,., +14^« -6^,.3 +^,.4!' - (38 

Coefficients set for multipliers 402-1 to 402-8 are same as 
coefficients of respective ^ presented by formula (38), and values of 
the coefficients reflect the corrections conducted for three times. 
5 Further, adding operation for ^ is conducted by the adder 24, and the 
added value is squared in a square device 25. Here, provided that ni-4 
to ni+4 satisfy probability distribution of Gaussian distribution in 
formula (37 ) and average power is presented as (na) ^, it can be expressed 
by formula (39) : 
10 [formula 39] 

<=f— 1 SSSn/ (39 

Here, assuming that a = (1/128)^ (858), noise power ("N+I other") 
is presented as: 
[formula 40] 

As such, according to the present embodiment, the estimation 
accuracy of the noise power can further be achieved, by adding corrections 
conducted for three times during the process for estimating noise power. 

20 (Fourth Embodiment) 

While the third embodiment illustrates a case of conducting the 
correcting operations for three times during the process for estimating 
noisepower, in the configuration of using ^ that are obtained by aver aging 
^ of adjacent twosubcarriers when noise power of subcarrier is calculated. 
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the fourth embodiment will illustrate a case of conducting the correcting 
operations for three times during the process for estimating noise power 
in the case of using ^ that are obtained by averaging ^ of adjacent 
four subcarriers • 

5 Fig. 9 is block diagram showing an internal configuration of 

a noise power estimation section according to the fourth embodiment 
of the present invention. In this figure, same numeral number as in 
Fig. 5 is assigned to an element that is commonly appeared in Fig. 5, 
and the detailed explanation thereof is omitted. The point that Fig. 

10 9 is different from Fig. 5, is that the number of the delay devices 
and the multipliers are respectively increased, .and coefficients for 
multipliers 502-1 to 502-16 and a multiplier 503 are changed. 

Arithmetic operations provided by the noise power estimation 
section shown in Fig. 9 will be described using formula. 

15 Here, specific arithmetic operations will be described focusing 

on a subcarrier i. However, since the operations until formula (31) 
is derived in this example is the same as in the second embodiment, 
description will be made for operations after obtaining formula (31) . 
In the fourth embodiment, additional correction is further added to 

20 formula (31) , and thereafter, similarly as each of the embodiments stated 
above, expansion and approximation are carried out to obtain the 
following formula (41) : 
[formula 41] 



= 



25 -366«,^2 +174/7,^3 +62«,^4 -SSw,^^ -6«,^ H-6«,,7 ^n,^^ (41 

Actually, the calculation of this formula (41) is carried out 



by replacing n with E,. 



[formula 42] 



<=[-^] 1-^,^ -6^/-7 +6^,., +58^,., -62^,., -174^,.3 + 366^, -238^, 
+ 238^,,, -366^,,, +174^,,3 +62^,^ -58#.,, -6^,^ +6^,,, +^,,,f 
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Coefficients set for multipliers 502-1 to 502-16 are same as 
coefficients of respective E, presented by formula (42), and values of 
the coefficients reflect the corrections conducted for three times. 
Further, adding operation for ^ is conducted by the adder 24, and the 
added value is squared in a square device 25, 

Here, provided that ni-s to ni+s satisfy probability distribution 
of Gaussian distribution in formula (41) and average power is presented 
as (na)^, it can be expressed by formula (43): 
[formula 43] 



Here, assuming that a = (1/2048)^ (456314), noise power ("N+I 
other") is presented as: 
[formula 44] 



As such, according to the present embodiment, the estimation 
accuracy of the noise power can further be achieved, by conducting 
corrections for three times during the process for estimating noise , 
in the case of using ^ that are obtained by averaging ^ of adjacent 
four subcarriers when noise power of subcarrier is calculated. 




(43 




(44 
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(Fifth Embodiment) 

While the second embodiment illustrates a case of adjusting number 
of subcarrier employed in the averaging operation, or in other words, 
a case of adjusting number of pilots arranged along the frequency 
5 orientation, based on the measurement result of a delay spread, which 
is measured by estimating the condition of frequency selective fading, 
the fifth embodiment of the present invention will illustrate a case 
of adjusting number of pilots arranged along the time orientation 
employed in the averaging operation. 

10 Fig, 10 is a block diagram showing a configuration of a receiving 

system according to the fifth embodiment of the present invention. 
However, same numeral number as in Fig. 6 is assigned to an element 
that is commonly appeared in Fig. 10 and Fig. 6, and the detailed 
explanation thereof is omitted. The point that Fig. 10 is different 

15 fromFig. 6, isthatthedelayspreaddeterminationsection201isreplaced 
with a fD estimation section 1001 and that the noise power estimation 
section 202 is replaced with a noise power estimation section 1002. 

The f D estimation section 1001 estimates Doppler frequency based 
on signal (received signal) that is output from a radio processing section 

20 11, and the estimated Doppler frequency is output to the noise power 
estimation section 1002 . 

The noise power estimation section 1002 changes number of PICK 
arranged along the time orientation employed for the in-phase adding 
operation, on the basis of the estimated results obtained by the fD 

25 estimation section 1001. More specifically, since the drift of fading 
is gradual when Doppler frequency is lower, and thus the accuracy for 
estimating noise power can be improved by increasing number of PICH 
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arranged along the time orientation employed in the in-phase adding 
operation. On the contrary, since the drift of fading is intense when 
Doppler frequency is higher, an increase of the error generated by the 
difference in channel variation can be prevented by decreasing number 
5 of PICH arranged along the time orientation employed in the in-phase 
adding operation. 

According to the present embodiment, number of PICH arranged 
along the time orientation employed in the in-phase adding operation 
is changed based on Doppler frequency, so that the influence to the 
10 fading is reduced, and the accuracy for estimating noise power is further 
improved . 

As have been described above, according to the present invention, 
the accuracy for estimating the noise power can be improved by repeatedly 
correcting for multiple times the differences in the channel variations 

15 in the different subcarriers when noise power is estimated, even if 
the correlation between adjacent subcarriers is low due to the frequency 
selective fading. 

Although each of the above preferred embodiments describes that 
corrected values of noise power per one subcarrier, which are obtained 

20 for the whole subcarrier wave, are cumulative-added, and the obtained 
values are averaged to estimate noise power, the present invention is 
not limited thereto, and corrected values of noise power per one 
subcarrier, which are obtained for a desired subcarrier, may also be 
cumulative-added, and the obtained values may be averaged by using the 

25 desired subcarrier. 

In addition, the multiplex system for pilot signal (PICH) 
described in each of the embodiments stated above may be any multiplex 
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system, such as time multiplex, f requencymultiplex, scatteredmultiplex, 
code multiplex and the like. In addition, the method for estimating 
noise power and the noise power estimation apparatus of the present 
invention may be applied to Multiple Input Multiple Output - Orthogonal 
5 Frequency Division Multiplexing (MIMO-OFDM) or transmission diversity . 
It is considered that, in the MIMO-OFDM or the transmission diversity, 
there are basically a lot of scatters, and the propagation paths for 
each of the antenna are influenced by highly-selective fading, and thus 
received power considerably fluctuates between subcarriers. 

10 Thus, the estimation accuracy of noise power can be improved 

by applying method for estimating noise power and the noise power 
estimation apparatus of the present invention. 

The present application is based on Japanese patent application 
No. 2003-038935, filed February 17, 2003, the whole contents of which 

15 are hereby incorporated by reference. 

Industrial Applicability 

The present invention is preferable to be applied to a radio 
20 receiver. 



